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O |' Abstract 

Prompted by the recent more precise determination of the basic cosmological parameters and growing 
evidence that the matter-energy content of the universe is now dominated by dark energy and dark matter 
1 we present the general solution of the equation that describes the evolution of density perturbations in 

the linear approximation. It turns out that as in the standard CDM model the density perturbations 
grow very slowly during the radiation dominated epoch and their amplitude increases by a factor of 
^ ' about 4000 in the matter and later dark energy dominated epoch of expansion of the universe. 

On 
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O ■ 1 Introduction 
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During the last decade the possibilities of observing the universe have substantially increased. Larger and 
better telescopes, more sensitive CCD cameras and high resolutions spectrographs have been built. These 
advances enabled more thorough investigation of the content of the universe and more accurate determination 
of the basic cosmological parameters. Recently undertaken project of measuring the Hubble constant and 
the deceleration parameter by using the type la supernovae as standard candles lead to surprising results. 
(2jT)! It turned out that the distant supernovae appear to be fainter than expected and in order to reconcile 
observations with theoretical expectations it is necessary to introduce the cosmological constant or dark 
energy as the dominant constituent of the present universe. The origin and nature of the dark energy is not 
know and therefore in what follows we assume that it can be associated with the cosmological constant. 

By now there is a solid evidence that large fraction of the total mass of galaxies and clusters of galaxies 
is in the form of non baryonic dark matter of unknown nature and composition. Existence of dark matter 
has been also confirmed by independent studies of the motion of galaxies and their spatial distribution and 
results of high redshift galaxy surveys (for example, 2dF and SDSS). 

Recently released data on the anisotropy of the cosmic microwave background radiation by the WMAP 
team when combined with observations of the type la supernovae at high redshifts and observed properties 
of clusters of galaxies and galaxy surveys are all consistently determining the basic cosmological parameters 
leading to what some already call the Standard Cosmological Model pQ. 

It is customary to express the density of different components of the universe in terms of the critical 
3H 2 

density g cr it = , where H is the present value of the Hubble constant and G is the gravitational 
o7rG 
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constant. So we have 



8nG 8nG kc 2 Ac 2 

n m = ^ Qm , r = ^ ei , k = -j^, and tt A = (1) 

where g m and g r are correspondingly the present density of matter and radiation including relativistic 
particles, k is the curvature parameter (+1, 0, —1), a — the present value of the scale factor, and A — the 
cosmological constant (dark energy). 

The fundamental parameters of the most popular cosmological models consistent with all sets of obser- 
vational data, can be taken as 

n m h 2 = 0.134 ± 0.006, n A = 0.73 ± 0.04, ^ k « 0, 

n r h 2 = 2.47 x 10" 5 , and H = 71 ± 4 km-s^-Mpc" 1 = lOOh km-s^-Mpc^ 1 . 

It is interesting to note that Ob/i~ 2 = 0.023 ± 0.001, and it is much less than O m and moreover only small 
fraction of baryons is luminous Oi um = (0.003 ± O.OOl)/^ 1 . It has been established that neutrinos are not 
massless but they have a small rest mass of the order of 0.02 eV and therefore if they decoupled from the 
state of thermodynamic equilibrium while relativistic their contribution to the density of relativistic particles 

is Vl^ti 2 = — — where rrn is the mass of z-th neutrino specie or 0„ ~ 6 x 10~ 4 . 
94 ev 

Prompted by these new results and in particular by the growing evidence that the cosmological constant 
(dark energy) is different from zero we want to analyze evolution of density perturbations in a model universe 
filled in with matter, radiation and the cosmological constant. 



2 Evolution of density perturbations 



We assume that the background universe is homogeneous and isotropic and its geometry is described by the 
Friedman-Robertson- Walker line element 



ds 2 =dt 2 -a 2 (t) 



dr 2 



1 — kr 2 



r 2 (d6 2 + sin 2 9dip 2 ) 



(2) 



where a(t) is the scale factor, k = +1, 0, —1, and we use units such that the velocity of light c = 1. The 
Einstein field equations with the hydrodynamical energy-momentum tensor describing the distribution of 
matter and the cosmological constant A reduce to: 



2 

a la' 



2-+ (-) = -8ttG P + A, (3) 

k 8ttG 1 



a \a 



- +- = ^7-2+o A - ( 4 ) 

a J a 2 3 3 

We assume that the universe is filled in with non relativistic matter (dust) and radiation (photons, 
neutrinos and other relativistic particles) but the non relativistic matter does not interact with radiation. 
We understand that this is a rather drastic assumption and it restricts the applicability of our model to 
the post decoupling period. The nature of dark matter particles — the dominant constituent of matter is 
not known. From laboratory and astronomical observations it only follows that they interact very weekly 
with baryons and photons and therefore dark matter particles decoupled from the state of thermodynamical 
equilibrium much earlier than baryons, quite possibly at z > 10 6 . 
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We can therefore decompose the total density of matter-energy into dust and radiation and so 

Q = Qr + £>m, (5) 

and 

P = 3^r- (6) 

In our model matter and radiation expand adiabatically and we have 

g r a 4 = const, (7) 

and 

QmO 3 = const. (8) 

The Friedman equation <@J supplemented by Q, JJJ and J^J form a closed system of equations. 
Unfortunately when g r , g m , and A are all different from zero the Friedman equation can not be solved 

analytically. However using the standard notation H = — and introducing the redshift parameter 1 + z = — 

a a 

equation QJ can be transformed into 

h 2 (z) = Hi (n,(i + zf + n m (i + z f + n k (i + zf + n A ) . (9) 

It is now apparent that if f2 r ^ then sufficiently early as z — > oo the density of radiation determines the 
expansion rate of the universe. 

The problem of evolution of small density perturbations in the Friedman universe was thoroughly dis- 
cussed by Lifshitz in 1946. The general relativistic equations describing evolution of small perturbations 
have been extensively studied since then. Here we are interested in evolution of density perturbations in the 
pressureless (dust) component of matter. Hypersurfaces orthogonal to the dust matter flow coincide with 
the equal proper time hypersurfaces, therefore the perturbation equation has the same form in both, the 
Newtonian theory and relativistic gauge- invariant formalisms El HI HI HI 

A + 2HA - AttGq^A = , (10) 

S Q 

where by capital A we denote the density contrast A = — measured on the flow orthogonal hypersurfaces, 

Q 

and dot stands for the derivative with respect to time t. Since the time dependence of H and g m is in the 

general case not known, instead of time we will use another independent variable x = — — - so x is changing 

a 

dx 

from to 1. Using the relation — = xH and J2J we transform the equation (| 1 Oil into 

3 3 

x(n r + n m x + n k x 2 + n A x 4 )A" + (o r + -n m x + 2n k x 2 + 3n A x 4 )A' - -n m A = o, (11) 

where the prime denotes differentiation with respect to x. 

The differential equation l(TT|l has a very special form, the coefficients of the second and first derivative 
are polynomials of the independent variable x. In turn the coefficients of these polynomials form a vector 
Q = {il m , 17 r , ilk, ^a} in the 4-dimensional parameter space. The Friedman equation <@J constrains these 
parameters imposing the condition that 

n r + n m + Q k + Q A = 1. (12) 
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In the parameter space this relation defines a hyperplane. It means that only three of the parameters are 
independent. The polynomial multiplying the second derivative in the equation (1 1 XII determines the number 
and positions of singular points of this equation in the finite domain. In the generic case this equation has 
six regular singular points including one at infinity. For certain values of fis some of the singular points 
became complex. To study the general properties of this equation it is therefore natural to consider the 
independent variable as complex. In the complex domain the equation 11H is of Fuchs type with six 
regular singular points. General solutions of such differential equations are practically unknown. We will 
seek and investigate solutions of this equation in two cases: when fl\ = equation (11 111 reduces to the 
Heune equation (HI El EDI E] with four regular singular points, and in the general case when fi A ^ the 
equation l(TT|l has six regular singular points. 

3 Solutions with £7a = 

When Oa = the equation ifTTJl assumes the form 

3 3 
x(fl T + n m x + tt k x 2 )A" + (fi r + -tt m x + 2tt k x 2 )A' ~ ~fi m A = . (13) 



-fi m T - 4fi k fi r 

This equation has four regular singular points at x\ = 0, X2.3 = and X4 = 00. 

General discussion of solutions of this equation in the parameter space {fi m , fi r , fik} will be provided after 
we present some of the particular solutions. 

• The case fi — {fi m , 0,0,0}. Exact solutions of equation ijlSj) are known for particular values of the 
parameters fi — {fi m , fi r , fik, ^a}- When fi — {fi m , 0,0,0}, the equation lllll can be easily solved and we 
get 

A = ciA + + c 2 A_ = c\x + c 2 x 5, (14) 

where by c\,C2 we denote arbitrary constants not only here but throughout this paper. It is apparent that 
A + (x) does vanish at x = and increases with increasing x. It describes the growing mode of density 
perturbations. The other solution A-(x) is singular at x — and decreases with increasing x. It describes 
the decaying mode of density perturbations. This solution was know already to Lifshitz in 1946 1121 . 

• The case £1 = {O m , O r , 0, 0}. Slightly more general solution of the equation 111 -ill is obtained when 

fi = {n m ,fir,o,o} ESQa 



A = n A_ + ' -.A. = r, I fi r + ^fi m x 



c 2 - 



15n ™ f 3Vfi r + fi m * - -1= ( fi r + |fi m aA In V " r + ) . (15) 



4 n 2 V V 2 / Vfir + n m x + VW,- 



As in the previous case A + (x) describes the growing mode of density perturbations, but now A + (0) = 
fi r . The other independent solution A_(x) is singular at x = and decreases with increasing x. The 
normalization of these solutions has been chosen in such a way that when fi r — * 0, A + and A_ tend to the 
corresponding modes of the solution ifTHl . 

• The case fi = {fi m , 0, fit, 0}. Another exact solution is obtained when the cosmological constant is zero 
and radiation is not present, in this case the non zero parameters are fi = {fi m , 0, fi k , 0} and the solution 
can be written in the form [l5llfi| 

A = ciA+ +c 2 A_ 
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, j\ I +^L X x-i. (16) 

As before A + (x) describes the growing mode of density perturbations, in this case A + (0) = and asymp- 

5 f2 

totically for large x it tends to lim A+(x) = - 7^. The second independent solution A_(x) is singular 

at x = and it decreases to zero when x — > 00. When fit — > then both modes of density perturbations 
A + and A_ tend to the corresponding modes of the solution ifHjl . 

• The case fi = {2^/^^■(l — y/Sh), f2 r , f2k, 0}. Another particular solution of the evolution equation (THfl) in 
a form of first several terms of a series expansion has been given by Guyot and Zeldovich |Tj[j] . This solution 
corresponds to the following choice of the parameters ft — {2-^/^(1 — f2 r , ^k, 0} and for this choice of 

parameters two singular points of the equation 111 -it X2 = X3 overlap. In this case the equation 111 -it assumes 
the form 

x(2n r + n mX ) 2 A" + 2{n r + n m x)(2Ci T + n mX )A' - 6n T n m A = o. (17) 

This equation can be solved exactly and the exact solution is expressed by hypergeometric functions 
A = d A+ + c 2 A_ 



ci f2 r 2 Fi [ Vs i, — V3 i, 1, ■ 



+ C 2 



2 F 1 [l + V3i,l-VSt,2, 



sinh[\/37r] 2fi r + fl m x \ ' ' 2f2 r + £l m x 



(18) 



Despite the appearance of imaginary unit the solutions are real because the first two arguments in the 
function 2 F\ are complex conjugate of each other. A + describes the growing mode with A + (0) = f2 r and 

asymptotically for large x we get lim A + (x) — fi r — — ^= . The second independent solution represents 

the decaying mode, it is singular at x = and when x grows to infinity this solution tends to zero. 

• The case fl = {ri m , fl r , 0}. Let us now return to the general discussion. From the expressions 
determining the positions of the regular singular points x 2 and xz it is apparent that in this case the 
parameter space is reduced to the plane of {fi r ,fi m }. Physically interesting is only a section of this plane 
described by f2 r € [0, 1] and Q m G (0, 1] shown on Fig.Q] On this figure there are three regions: A, B and 
C . The boundary curve between the regions B and C represents the condition that two singular points 
x 2 and X3 overlap. This curve is described by the equation Qj^ — 4ilk£l r = 0. This equation combined 
with the constrain 11211 with fl\ = leads to an explicit equation of this curve Cl m = 2y/Q^(l — V^r)- 
Therefore the solution presented by Guyot and Zeldovich Jgj is very special and in the parameter space 
{Q r , fi m } it is represented by a curve. The constrain l(T2|l with the additional assumption that Qa = ^k = 
reduces to an equation of a straight line Vt m + Cl x = 1, which is a boundary between the regions A and 
B. Therefore the regions B and C represent models with ilk > and the region A represents models with 
fik < 0. In what follows we will concentrate on the region B (shadow region), as it is the most interesting 
from the point of view of the present day observational data. The general equation describing evolution 
of the density perturbations l(T3j) in the region B with Q = {fi m , f2 r , Qk, 0} can be turned into the Heune 

(Q m + JUl. - 4n r n k )x 

equation by the following transformation of the independent variable £ = — , 

2fi r + (fi m + y/O^ - 4Q T Ct k )x 

which shifts the positions of singular points of equation 11-311 form {x±,X2, X3, X4} to {0, d, 00, 1}, where d is 
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Figure 1: The {f2 r ,f2 m } plane with various domains and singularity points, 
defined below. Then the parameters of the Heune equation are given by: a = |, @ = 0, 7 = 1, S = 0, e = ± 
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a, 



> 1, 9 



fin 



and the equation assumes the canonical form of the 



Heune equation IIA.ljl . Solutions of this equation near the singular point £ = can be written in a form of 
series l!A.4|l and l|A.8j) 



A = ci A + + c 2 A_ = ciF(<2, g; a, /3, 7, <5; ^) + c 2 G(d, q; a, (3, 7, 5; f ) , 



(19) 



where F and G represent two independent solutions given in Appendix A. As in the previous cases A + is 
regular at £ = and represents the growing mode and A_ is singular at £ = and describes the decaying 
mode. Solutions of the equation (17) with the redshift parameter z as an independent variable in a form of 
a series expansion have been given by Rozgacheva and Sunyaev [17] but they do not discuss convergence of 
these series. 



4 The general case when £7a 7^ 

• The case Q = {O m , 0, 0, Let us first consider some special cases of particular values of the parameters 
such that exact solutions of the evolution equation of density perturbations 11 111 are known. We start with 
a case when the universe is flat and without radiation, so fl — {f2 m , 0, 0, ^a}- In this case the general solution 
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of the equation 111 111 can be written down in the form 



1, 11 A 3 \ , L , A 3 _i 



5 •) o.l 1 / 3 

A + (x) is monotonically growing and asymptotically it tends to lim A + (x) = T 

6y 7f 



2 




5 






r 






3 




6 





A = Cl A + + c 2 A_ = c lX 2Fl 1, T , - j + c 2 ^ 1 + (20) 

(see P-9J). A + (x) describes the growing mode of density perturbations, in this case we have A + (0) = and 

"0_ A " 

_o n 

The second independent solution A_(x) is singular at x = and it decreases with increasing x asymptotically 
tending to lim A_(x) = When Oa — > the solutions A + and A_ tend to 11411 . The solution 

x— too V O m 

describing the growing mode has been found by Eisenstein (2DJ and written in the form of an elliptic integral, 
it turns out that it can be written in terms of hypergeometric function what is more convenient to handle. 
In j2] the solution A + (ir) is written in terms of the incomplete beta-function. 

• The case = {O m , 0, Ok, Oa}. When radiation is not present the parameter space is spanned by 
O = {O m , 0, Ok, Oa}- In this case the general solution of equation ifTT)) can be written down in the form JSj 

A = ci A + + c 2 A_ = Cl] /1 + -x + -x* * J {nm + QkX + QAx3)3/2 d X 

+ c 2 \ 1 + —x + —x 3 x 2. (21) 

The integral appearing in i'2l\) can be expressed in terms of the Jacobi elliptic functions and its explicit form 
is given in Appendix C. 

• The case = {O m , O r , Ok, Oa}. The most general case, when O = {O m , O r , Ok, Oa}, includes many 
cosmological models (see Appendix D) and therefore in what follows we will limit our discussion to the case 
when Oa > 0. From the weak energy condition for matter and radiation it follows that O m + O r > what 
with Oa > and Ok > restricts the parameter space to a cube with a size length equal to one [2] 

0<O m + O r <l, 0<O k <l, 0<O A <l. (22) 

This cube then represents only models of the type (H) in Fig. (see also Appendix D). In this general case 
the evolution equation 11 111 is of Fuchs type with six singular points. Positions of the proper singular points 
are determined by the polynomial in front of A"(x), so we have 

x(O r + O m x + O k x 2 + A a; 4 ) = 0. (23) 

From the imposed conditions Il22t it follows that the equation 112311 does not posses positive roots. It is 
apparent that one root is x% — 0. The other roots are either all complex or two are complex and two are 
negative. Using the standard technique of solving 4-th order algebraic equations |2S| we find that 



Xi 


= 0, 




X2 


1, 

= ^("i +vi + 


Wl), 


X3 


1, 

= 2^ Ul ~ Vl 


+ Wl) 


X4 


1, 


Wl), 


a?5 


1, 


- Wl) 


x 6 


= oo, 
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where we have added the singular point at infinity and 

1/2 



Ml 



Va 2 + b 2 + a 



Va 2 + b 2 



1 1/2 



Va 2 + b 2 + a 



1/2 



•\/a 2 + & 2 — a 



1/2 



wi = — \ — 2a — 2 



Ok 



b = - 



1/3 



-q/2 + Vv + -q/2-Vv 

1/3 



1/3 
1/3 



2 ft k 



while D = q 2 /4 + p 3 /27, p 



-q/2 + VV q/2-VV 

fll + 12fi r fi A 20g + 27n^fi A - 72ft r ft k Sl A 



27fi^ 



We ordered the roots in 



such a way that x 2 is the closest one to x\. With such order of the roots the homographic transformation 



x - x 2 



moves the singular points {x\, x 2 , x 3 , X4, 2:5, xq} of the equation lilt into {0, 00, di, d 2 , d 3 , 1}, 



3'3 



x± 



- , and c?3 



X5 



-. This transformation maps the interval x € (0, 00) into 



where g?i = 

X 3 - X 2 X4~ x 2 X5 — x 2 

£ G (0, 1). Under this transformation the equation 111 111 assumes its canonical form (IB. 11 with the following 
values of the parameters a = 5, (3 = 0, 7 = 1, 5 = —1, e\ = e 2 = e 3 = \. The additional parameters 

are given by q = -q , qi = -2q , q 2 = q where q = — — — — — — 2 • Tne general solution of the 

8S2 r + 6U m x 2 + ^UkX 2 

equation ltTT|l near the singular point £ = can be expressed in terms of the series l|B.4|) and IjB.QJl (see 
Appendix B) and we have 

A = ciA + + c 2 A_ = ciF(dx, d 2l d 3 , q, qx, q 2 ; a, f3,j, 6, ei, e 2 ; + c 2 G(dx, d 2 , d 3 , q, qi, q%\ a, (3, 1, 6, ei, e 2 ; £). 

(24) 

The solution A + represents the growing mode and A+(0) 7^ 0. The solution A_ represents the decreasing 
mode which is singular at £ = 0. 



Conclusions 

The general solutions of the equation that governs evolution of density perturbations derived and discussed 
above can be used to estimate the amplification factor of their amplitude. Unfortunately very little is known 
about properties of the dark matter particles and their interactions with the baryonic matter. We will assume 
that dark matter particles decoupled from the state of thermodynamical equilibrium at z w 10 s . This is 
quite conservative estimate but our general conclusions only weakly depend on this assumption. From the 
present values of radiation energy density and matter density f2 r w 5T0~ 5 and Q m w 0.7 it follows that until 
about z as 10 4 the universe was radiation dominated and from that moment to about z as 0.4 it was matter 
dominated and only recently its expansion rate became dark energy dominated. Using eq. J24II we calculated 
the amplification factor of the amplitude of density perturbations in a universe with O m = 0.3, f2 A = 0.7, 
n r = 5-10~ 5 and we got that A+(z = 10 4 )/A+(z = 10 8 ) 1.9, A + (z = 0)/A + (z = 10 4 ) « 3691 and 
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A + (z = 0)/A + (z = 10 8 ) 7012. Using eq. II24H we also calculated this amplification factor for an open CDM 
model with Q m = 0.3, Q r = 5-10" 5 and fi k = 0.7 and we obtained that A+(z = 10 4 )/A+(z = 10 s ) w 1.9, 
A+(z = 0)/A + (z = 10 4 ) « 2170 and A+(z = 0)/A+(z = 10 8 ) » 4122. This result confirms the well known 
fact (see that when the expansion rate of the universe is determined by the energy of radiation and 

relativistic particles the density perturbations grow very slowly (see eq. O)- The amplitude of density 
perturbations grows by a factor of about few times 10 3 in subsequent epochs. This puts strong constrains 
on the process of generation of the primordial density perturbations. The initial amplitude of the primordial 
density perturbations has to be sufficiently large to grow to A + w 1 at z w 10 in order to produce the 
observed large scale structure of matter distribution in the universe. This result only weakly depends on the 
existence of dark energy. It is now almost universally accepted that the primordial density perturbations 
have been created during the epoch of very early inflation. Unfortunately so far a sufficiently detailed model 
of inflation has not yet been proposed. Our results show that the constrains placed on the mechanism 
of generation of the primordial density perturbations by the observed large scale structure distribution of 
baryonic and dark matter only weakly depend on the existence of dark energy. 
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Appendix A 

A second order linear differential equation with four regular singular points z = {0, 1, d, oo} in its canonical 
form is known as the Heune equation [HI El El E], it has the form 

no + (2 + _i_ + m + m % ( «_ d) m - o, (a.i) 

where a, f3, 7, S, e, d and q are constants. The Heune equation is a natural generalization of the hypergeo- 
metric equation. The Heune equation and its solutions are determined by the Riemann P-symbol 

( 1 d 00 ~) 
P< a ;0 • (A. 2) 

(l-7l-Jl-e ^ J 

However this symbol does not fully determine the Heune equation and its solutions because the equation 
contains an additional parameter q. The parameters a, /3, 7, S, e are not independent and they are 
constrained by the Fuchs invariant 

l + a + (3 = 1 + & + e. (A.3) 

The solution of the equation i|A,lll , which is regular at £ = is known as the Heune function and it is 
defined by the series [lQ] 

00 

F{d, q; a, ft 7, *; = 1 + E Zu^C, (A.4) 

n=l "-wJnO 

where (j) n = 7(7 + 1) • • • (7+n— 1) is the Pochhammera symbol defined for 7 7^ 0, —1, — 2 . . . The coefficients 
c„ of the series l|A.4j) are determined by a recurrence relation of the third order 

Cn = {(n-l){a + fi- S + n-l + (7 + 5 + n- 2)d] + q} Cn-i- (o ! +n-2)(/3+n-2)(7+n-2)(n-l)dc n _2, 

(A.5) 
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for n > 2, with the initial condition 

c = 1, ci = q. (A. 6) 



The convergence radius of this series l|A.4|l is determined by the distance from the point £ = to the nearest 
singular point (the radius of convergence R = 1 for \d\ > 1 or R = \d\ for \d\ < 1). The second independent 
solution near the point £ = with the exponent 1 — 7 (7 ^ 1) has the form 

e~ 7 F(d, Ql ; 1 + a - 7, 1 + (3 - 7, 2 - 7, £; £), (A.7) 

where qi = g + (1 — 7)(1 + a + (3 — 7 + (d — 1)5). When 7=1 both solutions IIA.4J1 and l|A.7j) are identical 
and then the second independent solution contains a logarithm (see ^l]) and it has the form 

00 

G(d, q; a, (3, 1, 6; = F{d, q; a, (3, 1, 5; ln£ + ]T (A.8) 

n=0 

The coefficients 6„ of the series l|A.8j) are determined by the recurrence relation 

b n = -^{[q+(l + a + (3-S+(S + l)d+{d+l)(n-2)){n-l)]b n - 1 
- [a/3 + (a + (3){n - 2) + (n- 2) 2 ] 6„_ 2 - 2^^c n 
+ [l + a + (3-5+(S + l)d-(d + l)(l-2(n- ^l ^-i^. 1)|]a ^-i 
-[ Qi + /3 + 2(n-2)] d „_ 2[(n 1 _ 2)!]2 c„_ 2 }. 

for n > 2, with the initial conditions 



(A.9) 



6 = 0, b 1 = ~[a + (3 + S(d-l)-2c 1 ] . (A.10) 

Appendix B 

A second order linear differential equation with six regular singular points z = {0, 1, d%, d%, d$, 00} can be 
written down in the following canonical form 

+ ^ + + A + A) ™ + ^Xt-%U^ - 

and it is a generalization of the Heune equation. With the equation IIB.lj) one can associate the Riemann 
P-symbol 

{0 1 di d% d3 00 I 
a ;£ \ . (B.2) 
1 - 7 1 - 8 1 - ei 1 - £2 1 - e 3 J 

Because the equation l|B.lj) contains three additional parameters {q, gi, 52} therefore the Riemann P-symbol 
IIB.2jl does not uniquely determine the form of the differential equation. The Fuchs invariant puts additional 
constrains on the parameters that appear in the equation a, j3, 7, S, ex, €2, £3 imposing the relation 

l + a + /3 = "( + 5 + e 1 +e 2 + e 3 . (B.3) 
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The solution of the equation (IB . If) which is regular at £ = is a function that is a generalization of the 
series (IA.4I . it is given by 

oo 

F(d u d 2 , d 3 ,q, q u q 2 ;a, 0,j, 5, e x , e 2 ; = 1 + n \^)^d^d^ ' (B ' 4) 
The coefficients c„ in the series HB.4|) are determined by a recurrence relation of the fifth order 

C n = P n Cn-l + QnCn-2 + R n C n -3 + S n C n -4, (B-5) 

where 

P n = -q+(n~ l)(<2 2 d 3 ( 7 + e\ + n - 2) + di(d 3 ( 7 + e 2 + n - 2) 

+ tf 2 ( 7 + n - 2 + d 3 ( 7 + S + n - 2) + e 3 ))), 
Q n = -d!d 2 d 3 {n - 1)(7 + n - 2){q x + d 3 {n - 2)( 7 + e x + e 2 + n - 3) + d 2 {n - 2)( 7 + e x + n - 3 

+ d 3 ( 7 + S + e x ) + e 3 + n - 3) + di(n - 2)( 7 + e 2 + n - 3 + d 3 ( 7 + 5 + e 2 + n - 3) + e 3 

+ d 2 ( 7 + <S + e 3 + n-3))), 
i?„ = d\dldl{n - 2){n - 1)( 7 + n - 3)( 7 + n - 2)(12d 3 + (1 + d 3 )(n 2 + ( 7 + ei + e 2 - 7)n) - jq 2 

- 3d 3 ( 7 + S + ei + e 2 ) + ne 3 + d 3 £n + rf 2 (n - 3)( 7 + S + a + e 3 + n - 4) 

+ di(n - 3)( 7 + £ + e 2 + e 3 +n - 4) - 3( 7 + ei + e 2 + e 3 - 4)), 
S n = - 3)(n - 2)(n- l)( 7 + n - 4)( 7 + n - 3)( 7 + n- 2)(a/3 - 4( 7 + 5 + £l + e 2 + e 3 - 5) 

+ n( 7 + (5 + ei +e 2 + e 3 +n - 2)), (B.6) 

for n > 4, with initial conditions 

co = 1, 
Cl = q, 

c 2 = q 2 + [y(d 2 d 3 + di(d 2 + d 3 + d 2 d 3 )) + Sdid 2 d 3 + d 2 d 3 e x + d x d 3 e 2 + d\d 2 e 3 ] q - ^ d\d 2 d 3 q x , 
c 3 = q 3 + [(2 + 3 7 )(rfid 2 + d\d 3 + d 2 d 3 + d\d 2 d 3 ) + 3did 2 d 3 + 3(d 2 d 3 ei + d\d 3 e 2 + did 2 e 3 )]q 2 

+ [2dl<if(7 + + 7 + ei) + 2d 2 (d 2 ( 7 + e 2 )(l + 7 + e 2 ) + + ^) 7 2 + 2e 2 (d 3 <5 + e 3 ) 

+ 7(1 + ( 2 + 4(1 + <5 + £2) + £3)) + d 2 ((l + d 3 + d 2 3 )-f 2 + djS(l + 5) + e 3 + 2<2 3 5e 3 + ej 
+ 7(1 + 2e 3 + dg(l + d 3 + 5 + 2d 3 <S + e 3 )))) + d 1 d 2 d 3 (~{2 + 37)91 + 2(d 37 (l + 7 + ei) 
+ d 3 ( 7 + 2ei)e 2 + cf 2 ((l + d 3 ) 7 2 + 2eid 3 <5 + e 3 ) + 7 (1 + e x + d 3 (l + <5 + ei) + e 3 ))))]g 

- 2 7 [((5 + l)d 2 d2d 3 + did 2 c^ + d 2 d 2 d 3 + d 2 ^^ + ^1^2^361 + d\d 2 d\e 2 + djdld^q! + d\dldlq 2 ] 

- 2-y 2 [(d\dld 3 + d\d 2 d\ + d x d\dl + d\d\d 2 z )q x + d 2 d\d 2 3 q 2 ]. 

(B.7) 

The convergence radius of the series llB.4ll is determined by the distance between the point £ = and 
the nearest singular point (the radius of convergence R = min{l, \di\, \d 2 \, |d 3 |}). The second independent 
solution near the point £ = with the exponent 1 — 7 (7 ^ 1) has the form 

C 1_7 F(di, d 2 , d 3 , q, q u q 2 ; 1 + a - 7, 1 + - 7, 2 - 7, S, e x ,e 2 ; £), (B.8) 

where q = q+{l-j)(d 2 d 3 ei + di(d 2 d 3 5 + d 3 e 2 + d 2 e 3 )), q\ = q x + (1 ~7)(d 3 (ei +e 2 ) + rf 2 (ei +d 3 (5 + ei) + e 3 ) + 
di{e 2 +d 3 (5+e 2 )+e 3 +d 2 (S+e 3 ))), q 2 = <7 2 -(l-7)(d 3 <5+ei+d 3 ei+e 2 +d 3 e 2 +e 3 +d 2 (^+ei+e 3 )+di(J+e 2 +e 3 )). 
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When 7 = 1 both solutions llB.4ll and IIB.7II are identical and then the second independent solution contains 
a logarithmic term and is given by 

oo 

G(di,d2,d 3 ,q,qi,q2;a,/3, l,S,e 1 ,e 2 ;0 = F(d 1 ,d 2 ,d 3 ,q,q 1 ,q 2 ;a,f3, l,8,e 1 ,e 2 ;0 ln£ + b n C- (B.9) 

n=Q 

The coefficients b n of the series l|B.8l) are determined by the following recurrence relation 

bn = PnK-l + QnK-2 + J^fin-S + S%- 4 + P^C n + Q^C„_i + R C n C n - 2 + S£c„_ 3 + ^C„_ 4 , (B.10) 

where 

p n = oi\ , [q + d 2 d 3 (n - l)(n - 1 + e x ) + d x {n - l){d 3 {n - le 2 )d 2 {n —1 + 5) + e 3 ))], 
n s d\d 2 d 3 

Qn = 2 , . . [gi + di(n - 2)(n - 2 + d 2 (n - 2 + * + e 3 ) + d 3 (n - 2 + 5 + e 2 ) + e 2 + e 3 ) 



+ d 2 [n - 2)(n - 2 + d 3 (n - 2 + 5 + ei) + ei + e 3 ) + d 3 (n - 2)(n - 2 + ei + e 2 )], 
a , [?2 + (n - 3)(di (1 + 5 + e 2 + e 3 
(n - 4)(n - 3)(1 + rfi + rf 2 + efe) + £3)], 



i$ = , [ga + (n - 3)(di(l + 5 + e 2 + 63) + d 2 (l + 5 + + e 3 ) + d 3 (l + 8 + e 1 +e 2 ) 

n z did 2 d 3 



S n = *T-j . [a(3 + (n - A)(n - 4 + 5 + ex + e 2 + e 3 )], 

n z d\d 2 d 3 



n(n!) 2 (did 2 d 3 )™' 

Q„ = -77 77. ^» [ rf i( rf 3(2n - 2 + e 2 ) + d 2 (2n - 2 + d 3 (2n - 2 + 5) + e 3 )) + d 2 d 3 (2n - 2 + El )], 

[(n - l)!nj 2 (dici 2 £i 3 )™ 

^ = -77 0M 12 !, , , , n _ 1 [di(2n-4 + e 2 + £3 + d 3 (2n-4 + g + e 2 ) + d 2 (2n-4 + J + 63)) 

[{n - 2y.n\ 2 (d 1 d 2 d 3 ) n 1 

+ d 2 (2n - 4 + ei +e 3 + d 3 (2n -4 + 5 + e ± )) + d 3 (2n - 4 + e x + e 2 )], 
SI = [{n 3 y n]2{did2d3) ^ 2 [di(2n - 6 + J + Ca + e3 ) + d 2 (2n - 6 + 5+ e x + e 3 ) 
+ d 3 (2n - 6 + 5 + ei + e 2 ) + 2n - 6 + e x + e 2 + e 3 ], 

1 —[2n -8 + 5 + £l +e 2 + e 3 ], (B.ll) 



[(n-4)!n] 2 (d 1 d 2 d 3 )™- 3 
for n > 4, with the initial conditions 
bo = 0, 

6 1 = , , , [d 2 d 3 ei + di(d 3 e 2 + d 2 (d 3 5 + e 3 )) - 2c x ] , 

«lCl2Ct3 

62 = \d 2 d 2 d 2 ^ 2d3ei + dl ( d 3 £2 + rf 2(d 3 <5 + £3))]? + d 2 d 3 ei(d 2 d 3 (l + ei) + 2did 3 e 2 + 2did 2 (d 3 5 + e 3 )) 

+d 2 (d 2 e 2 (l + £2 ) + 2d 2 d 3 e 2 (d 3 5 + e 3 ) + d^t^l + 5) + (1 + 2d 3 <5 + e 3 )e 3 )) 
-~[2q + d 2 d 3 ex + d 1 (d 2 d 3 5 + d 3 e 2 + d 2 e 3 )]ci - c 2 } , 

63 = 1Q8d l d 3 d 3 {3[dadsei + d x (d 3 e 2 + d 2 (d 3 <5 + e 3 ))]g 2 + 3[-d 2 d 2 ei(5 + 3e x ) - 2d 1 d 2 d 3 {d 3 (2e 1 + 2e 2 + 3eie 2 ) 
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+ d 2 (2d 3 S + 2ei + 2d 3 e l + 3d 3 <5ei + 2e 3 + 2 ei e 3 )) - d?(d^e 2 (5 + 3e 2 ) + 2d 2 d 3 (2e 2 + d 3 (2S + 2e 2 + 3<5e 2 ) 
+ 2e 3 + 3e 2 e 3 ) + d|(^<5(5 + 36) + e 3 (5 + 3e 3 ) + d 3 (A6 + 4e 3 + 6Se 3 )))]q 

- 12did 2 d 3 [d 2 d 3 ei + dx{d 2 e 2 + d 2 {d 3 5 + e 3 ))]gi + 6(d|d|ei(2 + 3ei + e?) + 3^d 2 d|eie 2 (l + e 2 ) 
+ 6df^4eie 2 (d 3 5 + e 3 ) + 3d?dSSd 3 ei(dt<J(l + (5) + 2d 3< 5e 3 + e 3 (l + e 3 )) 

+ 3di^^ei(l + ei)(d 3 e 2 + d 2 (d 3 5 + e 3 )) + d? (d|e 2 (2 + 3e 2 + 4) + 3d 2 d 3 e 2 (l + £ 2 )(d3<S + £3) 

+ 3d|d 3 e 2 (d|(5(l + 8) + 2d 3 Se 3 + e 3 (l + e 3 )) + d|(d|(5(2 + 3<S + S 2 ) + 3d 3 <5(l + <S)e 3 + 3d 3 <Se 3 (l + e 3 ) 

+ e3 (2 + 3 £3 + e 2)))) 

+ [-6q 2 + 3(d 2 d 3 (8 + 5ei) + did 3 (8 + 5e 2 ) + d x d 2 (8 + 8d 3 + 5d 3 S + 5e 3 ))q 

+ 3(-2d|d|ei(2 + ei) - 4d 1 d 2 d 3 (-2q 1 + d 3 eie 2 ) - 4did| ! d 3 ei(d 3 <$ + e 3 ) - 2d\{d 2 3 e 2 {2 + e 2 ) 

+ 2d 2 d 3 e 2 (d 3 £ + e 3 ) + d 2 (d 2 S(2 + S) + 2d 3 8e 3 + (2 + e 3 )e 3 ))))]ci 

- 3[q + d 2 d 3 ei + di(d 2 d 3 5 + d 3 e 2 + d 2 e 3 )]c 2 - 2c 3 } . 

(B.12) 



Appendix C 

To evaluate the integral 

x 3/2 

(n m + n kX + n A x3)V 2dx ' (C1) 

it is necessary to resolve the third order polynomial 

n m + n k x + n A x 3 = 0, (c.2) 

into elementary prime factors. The number of real roots of this polynomial IIC.2II is determined by the sign 
of its discriminant 

_ 4nj + 2m 2 n n A 

There are three possibilities: 
(a) V = 0. This occurs when il A — — 2 7n^ and then the polynomial IIC.2JI can be represented as 

n m + n k x + n A x 3 = — ^-(3n m - n k a:)(3n m + 2n k x) 2 = o. (c.4) 

' m 

In this case the integral l|C.l|l can be expressed by elementary functions and we get 
* (0) models (Fig.© for f2 k > and x € (0, 3fl m /fl k ) 



2n 2 ^m m - n iX {3n m + 2n k x) 2 2 nl /2 V 3^ m - n k x 



and 



V30nV^ (9^ n + 9Q k Q m x + 8Q 2 x 2 ) fl m I 37^ 

arctan \/^=; Fi — 1 (C.5) 



* (O) models (Fig.0 for O k < and s e (0, -§fi m /ft 



V3n mv /i (9Q 2 m + 9fl k n m x + 8^x 2 ) _ Q m arctanh / 3^ k x 
2^2 V3^ m - O k x (3fi m + 2fl k x) 2 202y^I\ arCan V 3f7 m -fi k a;' 1 ' j 
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(b) V > 0. This happens when ft A > and 4ft 3 + 27fi^fi A > or ft A < and 4ft 3 + 27S1^S1 A < 0. In this 
case the polynomial (C.2) has one real root and two complex ones. In this case the polynomial l|C.2Jl can 
be written as 



ft m + Q]<x + Qj\x 3 = Q A (x — a)(q + px + x 2 ) = 

1/3 



(C.7) 



v 1/3 


/ f2m 




) + 


Vv) 




V 20 A 



- a 



where a = (-|g- + v^ +{-^~V^) , q - 
is given by 

* (H) models (Fig.EJ for ft A > and 4ft 3 + 27fl A ^m > 0, and x e (0, oo 
2(-9q 2 a 2 + 3<?a(q - a 2 )x + 2(<? - a 2 ) 2 x 2 ) 
ft A /2 (-4g + a 2 )( 9 + 2a 2 ) 2 y/xjx - a){q + ax + x 2 ) 



and p = a. In this case the integral iC.ll 



18q 3 / 2 a 2 



yj (x — a) (q + ai + x 2 ) 



ft A /2 (4g - a 4 )(q + 2a 2 ) 2 yft(y/q{x ~ a 2 )y/q + 2a 2 x) 



18g 5 / 4 



a 



n\ ,2 {Aq-a i ){q + 2a 2 y/ i 



E (arccos(u), k) 



where u 



^Jq(x — OL)—x\J q+2a 2 



and k 



2q+a 2 



3q 3 / i a(-3 v r g+ y/g + 2a 2 ) 
n A /2 (4q-a 4 )(q + 2a 2 ) 7 / 4 

1/2 



F(arccos(w),fc) . (C.8) 



Ay/q(q+2a?) 



and by 

* (O) models (Fig.0 for ft A < and 4ft 3 ; + 27ft A ft 2 „ < 0, and x e (0, a) 

2((2< 7 + a 2 )a; + 3ga:) v /i 
(-ft A ) 3 / 2 (4<? - a 2 )(q + 2a 2 )y/(a - x){q + ax + x 2 ) 
18qa 2 



\/g + as + 



(-^ /2 )(4g - a 4 )(g + 2a 2 ) 3 / 2 y/a~^(y/q(a - x) + y/q + 2a 2 x) 



18g 5 / 4 c 



-ft A ) 3 / 2 (4 g -a 2 )(g + 2a 2 ) 7 /4 



E (arccos(u), k) 



where u = — ^=^= and k 

y/q(a - x) + Xyf q + 2a 2 



3g 3 / 4 a(3yg+ y/q + 2a 2 ) 
(~n A ) 3 / 2 (4q - a 2 )(q + 2a 2 ) 7 / 4 

9 1 !/ 2 
2g + a 2 



F (arccos(w), k) , 
(C.9) 



2 4 y /q(q + 2a 2 )_ 

(c) D < 0. This case occurs when ft A > and 4ft 3 + 27ft^ft A < or ft A < and 4ft 3 + 27ft 2 1 ft A > 0. Then 
the polynomial ljC.2l) can be written as 

ft m + n k x + n A x 3 = n A { x - a)(x - 0)(x - 7) = o. (c.10) 

When ft A < and 4ft 3 + 27ft 2 n ft A > we have 



ft^ 



-3ft A 



where <j> = arccos 20a , _ 3 q a 
When ft A > and 4ft 3 + 27ft 2 ft A < we have 







-3/2 



ft^ 



-3ft A 



cb IT 

3 + 3 



fti< 



-3ft A 



(J) IT 

3 ~ 3 



" =2W 3ft7 C0S 



<j) 7T 

3 ~ 3 



i = 2 ^3ft7 C0S 



3 + 3 



7 = -2 



-ftk 
3ft A 



■ cos 
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where 



20 A 



( - Q A 

\3U A J 



-3/2" 



* This case corresponds to (0) models (Fig. HJ| for Q A < and 4ft 3 . + 27ft 2 n ft A > 0. The roots of the 
polynomial IIC.lOll have been chosen so that a > > (3 > 7, and x € (0, a). In this case the integral 
IjC.ljl can be expressed by the elliptic integrals 

2(a 7 (a - 2/3 + 7) + {a{/3 - 2 7 ) + fa)x)y/x 



(a - /3)(/3 - 7) (a - 7) V(« - " ^ - 7)(-«a) 3 / 2 
27[(/3 + 7)(a 2 + ^7) + «(a 2 - 6/?7 + 7 2 )] 



(a _/3)( a - 7 )209- 7) 2 V Z 7(^ 
2/3 7 (-2a + /3 + 7) 



/3)(-ao 3 / 2 



i? (arcsin \A<, fc) 



(a - /3)(a - 7)09 - 7) V~7(a = P)(-^a) 3/ ' 



-F (arcsin y/u, k) 



(C.11) 



{a-f})x 
where zt = —7 -tt- and k = 

a{x — (3) 



-1 1/2 



(a - /3) 7 



* (B) models (Fig. EJ) for ft A > and 4ft 3 . + 27ft 2 n ft A < 0. The roots of the polynomial i((TTo}) have 
been ordered so that a > (3 > > 7, and x £ (a, 00). The integral lIC.lll is given in terms of the elliptic 
integrals 



2/3[(/3 + 7 )(q 2 + 0j) + a(f3 2 - 6/3 7 + 7 2 )] y/{x - a)(x - 7) 

[(a -/?)(«- 7 )(/3 -7)] 2 ^ /2 V*(*-P) 
2a/3 7 ((/3 + 7 )(a 2 + fry) + a(/3 2 - 6/3 7 + 7 2 )) 
[(a _ /3)( a _ 7 )09 - 7 )]2ft 3 / 2 x /a;(x-a)(x-/3)( a: -7) 
| 2[/3 2 7 2 (-2a + g + 7) + a 3 (/3 2 + 7 2 ) + « 2 Q3 + l){P 2 - Wl + 7 2 )]z 
[(a - /?)(a - 7 )0(3 - 7 )] 2 ft A / Vz(x - a)(a; - /3)(x - 7) 
4[ a 2 (/3 2 - /? 7 + 7 2 ) + /3 7 (/? 7 - qQ3 + 7 ))]x 2 
[( a _ (3)( a - 7 )(/3 - 7 )] 2 ^ A / V^ = ~ ~ l) 



2^[(f3 + 7) (a 2 + fry) + a(/3 2 - 6/3 7 + 7 2 )] 
(a - /3) 2 (a - 7 ) 2 (/3 - 7) 3 / 2 ft A /2 
2Va7(a - 2/3 + 7) 



£7 (arcsin */u, fc) 



(a _^ (a _ 7 )2 (/9 _ 7 )3/2 n 3/ 2 



.F (arcsin y/u, fc) 



(C.12) 



where u 



C3-7)(s- 
(a - 7)(x - /3) 



and fe 



g(g ~ 7) 
a{(3 - 7) 



1/2 



* (O) models (Fig.0J for ft A > and 4ft 3 . + 27ft 2 n ft A < 0. The roots of the polynomial (fCTTtTl have been 
ordered so that a > [3 > > 7 and a; € (0, (3). In this case the integral l|C.l|l can be also expressed by 
the elliptic integrals 

2a/3 7 [(/3 + 7) (a 2 + /3 7 ) + a(/3 2 - 6/3 7 + 7 2 )] 



[( a _ /3)( a _ 7 )(/3 _ 7 )]2r ] 3/2^ a , (a „ a;)(/ g„ x)(:E „ 7) 
2[/3 2 7 2 (-2a + /3 + 7) + a 3 Q3 2 + 7 2 ) + a 2 {P + l)ifl 2 - Wl + 7 2 )> 
[(a - /?)(a - 7)08 - 7 )] 2 ^a / V^(« " = ^ = 7) 
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4[g 2 (/? 2 - /fry + 7 2 ) + /fryQfry - g(g + 7 ))]x 2 
3/2 



[(a - /3)(a - 7 )(/3 - 7)]^r#^MM 



2a[a 2 (/3 + 7 )(a 2 + /? 7 ) + a(/3 2 - 6/? 7 + 7 2 )] y/(/3 - x)(x - 7) 
[(a - /J)(a - 7) 09 " 7)] 2 ^ A /2 " *) 

+ 2^ + 7)(^ + ^) + a(/? 2 -6/ 3 7 + 7 2 )] j; ^ ^ fc) 

(a-/3) 2 (a-7) 2 (/3-7) 3 / 2 r!A 
, 27(03 + 7 )(« 2 + /3 7 ) + ^ ~ 6/37 + 7 2 )]^ (arcgin ^ fc) 



V^(a-/3)(a-7) 2 (/3-7) 5 / 2 ^A /2 
2/3 7 (-2a + /? + 7) 



v^(a-/3)(a-7) 2 (/3-7)5/2^/ 2 



F(arcsinv^,/c), (C.13) 



where Hi = —. r, u 2 = —, r and k 

p{ot — x) p(x — 7) 

Appendix D 



P(a - 7) 



a(J3 - 7) 



1/2 



The standard classification of the Friedman cosmological models describing evolution of the universe that 
is filled in with dust, radiation and non zero cosmological constant is based on a technique used already 
by Friedman who studied the behavior of the curve x(t) = (the so called zero velocity curve). The zero 
velocity curve corresponds to the isoclines in higher dimensional dynamical systems. (However in our case 
the equation x(t) — does not determine positions of the fixed points. In the Friedman equation x(t) 
appears squared and therefore in the equation x(t) = F{x) the right hand side does not define diffeomorphic 
mapping at x = 0. 

In the classical papers ([211 EH, see also OI23EE1) the zero velocity curve is determined on the surface 
{x,Qa}- We use a more transparent classification on the surface {x, fl^} on which it is easy to identify 
solutions that correspond to phase trajectories in the classification of Belinsky and Khalatnikov [23|. From 
the Friedman equation © and © , © , © and the definitions of the Omega parameters © we get 

1 Pit) a , n m , ^ k , ^ 



H 2 x 2 {t) x A {t) x 3 {t) x 2 {t) 



where x(t) = In this case the equation that determines the isocline has the form 



ao 



-ft k = — !— + (l + x + x 2 )fl A - — ■ (D.2) 
x — 1 X 

The isocline has a vertical asymptote at x = 1. From the point of view of cosmological applications the 
most interesting is the region where S! m > and f2 r > 0. In Figs. EJEI and 21 the non physical region where 
x is complex is shaded and denoted as "forbidden". In Figs. El and 31 we represent respectively the allowed 
parameter space for A > 0, A = and A < restricted to the half plane representing universes with positive 
mass-energy density. 

In each of this figures the allowed cosmological models are represented by horizontal segments or half 
lines in the non shaded region that is bounded by the zero velocity curve. The models (O) are oscillating. 
For such models the universe starts its evolution form singularity and the universe expands until x = and 
then contracts to the final singularity. The models (B) are possible only for A > 0, they represent the non 
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-Ok 




forbidden 



(B) Lanczos 



(E) Einstein 



(E 2 ) 



(L) Lemaitre 



(P) parabolic (spatially flat) 



1 



(H) hyperbolic (negative curvature) 



de Sitter 

Figure 2: The {x, ilk} plane for fl\ > 0. 




(H) hyperbolic (negative curvature) 



Minkowski 

Figure 3: The {x, ^k} plane for fl\ = 0. 



singular Lanczos models. Such models describe the universe that starts expanding from a finite volume and 
expands to infinity. There is a possibility of bouncing when the universe initially contracts having a finite 
volume larger than the minimal possible for the given set of parameters. The half line that is tangent to the 
zero velocity curve represents two models Ei and E2 both approaching the static Einstein model (E) (the 
point of intersection) correspondingly in the forward and backward direction in time. The horizontal half line 
below the critical one represents the Lemaitre models. The x axis represents the flat models (parabolic (P) 
models in the Ehlers-Rindler classification [2Hj) and the half lines below represent universes with negative 
spacial curvature (hyperbolic (H) models in the Ehlers-Rindler classification [2Hj)- When A = the space of 
solutions contains only models (0), (P), (H) and the Minkowski space. When A < only oscillating models 
(0) are possible. 

Let us note that the zero velocity curve x = does not cross the line x — 0. The crossing condition 



Hax + il m x + fl Y = 



(D.3) 
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4 



3 



2 



1 



-1 




Figure 4: The {x, f^} plane for f^A < 0. 



does not have real solutions when all {fiA,O m ,O r } are positive. It means that the cosmological models 
with negative spacial curvature and positive cosmological constant independently on the values of the other 
cosmological parameters always start their evolution from the initial singularity and expand to infinity. We 
consider those models in our analysis of evolution of density perturbations. 
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